Scaling relations in the quasi-two-dimensional Heisenberg antiferromagnet 
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The large- A'^ expansion of the quasi-two-dimensional quantum nonlinear cr model is used in or- 
der to establish experimentally applicable universal scaling relations for the quasi-two-dimensional 
Heisenberg antiferromagnet. We show that, at A'^ = oo, the renormalized coordination number in- 
troduced by Yasuda et al. [Phys. Rev. Lett. 94, 217201 (2005)] is a universal number in the limit 
of J'/J 0. Moreover, similar scaling relations proposed by Hastings and Mudry [Phys. Rev. 
Lett. 96, 027215 (2006)] are derived at A'^ = cx3 for the three-dimensional static spin susceptibil- 
ity at the wave vector (7r,7r,0), as well as for the instantaneous structure factor at the same wave 
vector. We then use 1/N corrections to study the relation between interplane coupling, correlation 
length, and critical temperature, and show that the universal scaling relations lead to logarithmic 
corrections to previous mean-field results. 
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I. 



INTRODUCTION 



Many materials, such as copper oxides and other lay- 
ered perovskites, are known to be nearly two-dimensional 
magnets. While in certain intermediate temperature 
ranges these systems are well described by purely two- 
dimensional models, three dimensionality is restored at 
temperatures below an energy scale that is governed 
by the ratio between the interlayer coupling J' and 
the intraplanar exchange parameter J. For example, 
two-dimensional quantum Heisenberg antiferromagnets 
(AFs) do not support long-range coUinear magnetic or- 
der at any finite temperature according to the Mermin- 
Wagner theorem^ while real layered systems, such as 
La2Cu04, do. The anisotropy J'/J of a quasi-two- 
dimensional material can be determined from the spin- 
wave dispersion below the ordering temperature. It can 
also be determined from the measured ordering temper- 
ature provided one understands the two-dimensional to 
three-dimensional crossover that manifests itself in the 
dependence on J'/J of the ordering temperature. 

A common approximation for the ordering tempera- 
ture Tj, of a quasi-n-dimensional magnetic system is the 
random-phase approximation (RPA)ji^i^ which predicts 
that 
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(1.1) 



Here Xc"^ is the exact static susceptibility associated to 
the magnetic order for the underlying n-dimensional sub- 
system evaluated at the ordering temperature and z*-"-* 
is the coordination number of the n-dimcnsional subsys- 
tem. Yasuda et al. in Ref. have quantified the accuracy 
of the RPA in two steps. First, they computed the three- 
dimensional AF ordering temperature Tj^p with the help 
of a quantum (classical) Monte Carlo (MC) simulation of 
a spin-5 (5 = oo in the classical system) nearest-neighbor 
Heisenberg model on a cubic lattice with AF exchange 
coupling J' along the vertical axis and AF exchange cou- 
pling J within each layer of the cubic lattice. Second, 
they computed the two-dimensional static staggered sus- 
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ceptibility Xs evaluated at Tj^p from step 1 after switch- 
ing off J'. They thus showed that, for small J' / J , Tj^p 
is given by a modified random-phase approximation, in 
which the coordination number gets renormalized. 
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It turns out that the rcnormalization of the coordi- 
nation number z'^^ = 2 converges as J'/J ^ to a value 
that is independent of the spin quantum number S tak- 
ing values in 1/2, 1, . . . , oo. This fact motivated them to 
conjecture the universality of the effective coordination 
number C*^^) in the limit J'/ J — *■ 0. 

The results of Yasuda et al. were shown by Hastings 
and Mudry in Ref. H to reflect the so-called renormalized 
classical (RC) regime of the underlying two-dimensional 
subsystem. Hastings and Mudry predicted that if the 
two-dimensional subsystem is characterized by a quan- 
tum critical (QC) regime, then the effective coordina- 
tion number in the limit J'/J ^ is a universal func- 
tion of the ratio c/?^^/^'^' - a number of order 1 in the 
QC regime as opposed to a vanishing number in the RC 
regime - where c is the two-dimensional spin-wave ve- 
locity, Pj^p — 1/Tj^p, and ^'^^ is the two-dimensional 
correlation length. 

Hastings and Mudry also proposed universal scal- 
ing relations involving observables of the quasi-two- 
dimensional system only. One of these scaling relations 

is obtained from multiplying the static three-dimensional 
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spin susceptibility Xs evaluated at the wave vector 
(7r,7r, 0) and at the Neel temperature with J'. Another 
scaling relation can be derived by multiplying the in- 
stantaneous structure factor at the transition tempera- 
ture S*^^^ (tt, tt, 0; i = 0; T^p) with J' and the inverse Neel 
temperature l/Tj^p. This last universal relation has the 
advantage of being directly measurable with the help of 
inelastic neutron scattering. 

The aim of this work is to verify the universal scaling 
laws proposed above within the N = oo approximation 
to the quasi-two-dimensional quantum nonlinear a model 
(QNLSM). Our starting point is the two-dimensional 
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QNLSM (a two-space and one-time NLSM), which is be- 
Ueved to capture the physics of two-dimensional quan- 
tum Hciscnbcrg AF at low energies.— Assuming that 
the physics of a single layer can be approximated by 
the two-dimensional QNLSM, we introduce the inter- 
layer coupling following Refs. @, 0, M, and We 
then show that, in the N = oo approximation, the 
quantities J'x^^\t^f), J'^^^^tt, tt, 0; w = 0,T^p), and 
J'S*^^^ (7r,7r,0;i = 0',T^p)/T^p each converge to univer- 
sal scaling functions of the dimensionless ratio between 
the thermal dc Broglie wavelength of spin waves and a 
correlation length in the limit J'/ J 0. As in the strict 
two-dimensional limit, different regimes can be distin- 
guished depending on the value taken by this dimension- 
less ratio. The form of the nonuniversal corrections to 
the universal functions obtained in the limit J' /J 
strongly depends on these regimes. In the renormalized 
classical (RC) regime, which is dominated by classical 
thermal fluctuations, the nonuniversal corrections are rel- 
atively small for small J' /J, as is seen from MC simula- 
tions.— Near the critical coupling of the two-dimensional 
QNLSM, fluctuation physics is predominantly quantum. 
In this quantum regime, the universal constants obtained 
in the limit J'/J —^ depend on the way this limit is 
taken, i.e., on the value taken by the fixed ratio (14. 5p 
between the characteristic length scale of the quantum 
fluctuations and the characteristic length scale at which 
interplane interactions become important. When this ra- 
tio is a large number, the fluctuation regime is the quan- 
tum critical (QC) regime. When this ratio is of order 1, 
the fluctuation regime is the quantum disordered (QD) 
regime. When this ratio is much smaller than 1, the 
quasi-two-dimensional system remains disordered all the 
way down to zero temperature. 



The universal relations described above have recently 
been tested numerically by Yao and Sandvikiii They 
have performed quantum Monte-Carlo simulations of 
spin-1/2 quasi-two-dimensional systems in the RC, QC, 
and QD regimes. Their results qualitatively agree with 
the predictions of the large- expansion proposed here. 

We state the universal scaling laws in Sec. [TTl The 
quasi-two-dimensional QNLSM, which is the basis of our 
analysis, is defined in Sec. IIIII The N = oo approx- 
imation to the Nccl temperature is given in Sec. IIVI 
Scaling laws in the N ~ oo approximation are derived 
for the quasi-two-dimensional QNLSM in Sec. IVland for 
its classical limit, the quasi-two-dimensional NLSM, in 
Sec. IVII Conclusions are presented in Sec. IVIlIl while we 
defer to the appendixes for the derivation of the coun- 
terparts to these scaling laws in the case of the quasi- 
one-dimensional Ising model. Among the most impor- 
tant conclusions is an extension to finite A'^, where we re- 
late J', J, Tap, S}-'^'^^ and find scaling relations that differ 
from mean-field estimates by logarithmic corrections. 



II. UNIVERSAL SCALING RELATIONS 

We are going to formulate precisely the universal rela- 
tions discussed in Sec. HI It is instructive and necessary 
to first identify the relevant length scales of the system 
near the crossover between two dimensionality and three 
dimensionality. These will be used to construct the rel- 
evant dimensionless scaling variables whose values will 
allow us, in turn, to separate different scaling regimes. 

A. Characteristic length scales 

We start with the Hamiltonian describing a spatially 
anisotropic Hcisenberg AF on a cubic lattice, 

^ E ■ s^'. + ^'E E • (2-1) 

("') i i {jj') 

Planar coordinates of the cubic lattice sites are denoted 
by the letters i and i'. The letters j and j' label each 
plane. The angular bracket {ii') denotes planar nearest- 
neighbor sites. The angular bracket {jj') denotes two 
consecutive planes. A spin operator located at site i in 
the J th plane is denoted by j . It carries the spin quan- 
tum number S = 1/2, 1, . . . . The coupling constants are 
AF, J, J' > 0, and are strongly anisotropic, J ^ J' . 

We are interested in the thermodynamic properties of 
the system at the AF ordering temperature Tj^p. At this 
temperature, there exists a diverging length scale, the 
temperature-dependent correlation length 

e = 6%, (2.2) 

for AF correlations at the wave vector (7r,7r, tt). For a 
quasi-two-dimensional AF, the strong spatial anisotropy 
J ^ J' > allows us to identify the characteristic length 
scale 




(2.3a) 



where a is the lattice spacing of the cubic lattice and 
Z' is a multiplicative renormalizationJ^ In the N = oo 
approximation considered in this paper, = 1. 

We will omit this renormalization factor in the following 
until we treat 1/N corrections in Sec. IVII XI 
In a good quasi-two-dimensional Heisenberg AF 

< J' <C J =^ > a- (2.3b) 

It was proposed^iiS that ^' determines the crossover be- 
tween two and three dimensionality. Heuristically, the 
temperature dependence of the AF correlation length is 
nothing but that of the two-dimensional underlying sys- 
tem as long as ^^"^^ <€. C ■ Upon approaching T^^ 
from above, ^'^'^^ grows until .j'^' ~ . Below this temper- 
ature, three-dimensionality is effectively recovered and 
AF long-range order becomes possible. 
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In the regime of temperature above the ordering tem- 
perature T^p for which ^^^^ ^ a, the two-dimensional 
spin-wave velocity c is well defined and depends weakly 
on temperature. Together with the inverse Neel tempera- 
ture pj^p, we can then build the thermal de Brogiie wave- 
length cj3j^p. Equipped with the characteristic length 
scales c(3j^p and we can define three regimes, (i) The 
renormalized classical (RC) regime is defined by the con- 
dition 



(2.4) 



(ii) The quantum critical (QC) regime is defined by the 
condition 



a < c/3 



AF 



(2.5) 



(iii) The quantum disordered (QD) regime is defined by 
the condition 



(2.6) 



In the first regime (|2.4p . fluctuations are predomi- 
nantly two-dimensional and thermal. The planar cor- 
relation length diverges exponentially fast with decreas- 
ing temperature. This leads to a rather sharp crossover 
between two- and three-dimensionality. As a corollary, 
there will be small nonuniversal corrections to univer- 
sality for finite J'/ J. In the second, Eq. (|2.5p . and 
third, Eq. (|2.6[) . regimes, fluctuations arc predominantly 
two-dimensional and quantum. In the QC regime (|2.5p . 
the underlying two-dimensional fluctuations are quantum 
critical. In the QD regime (|2.6p . they are quantum dis- 
ordered. 

Common to all three regimes (|2.4p - (j2.6p is the fact that 
the lattice spacing is much smaller than the de Brogiie 
wavelength constructed from the two-dimensional spin- 
wave velocity and the three-dimensional ordering tem- 
perature. It is then reasonable to expect that planar 
fluctuations of the microscopic Hamiltonian (j2.ip can be 
captured by an effective low-energy and long-wavelength 
effective continuum theory. We choose this effective field 
theory to be the two-dimensional QNLSM.^^iii^ To ac- 
count for intcrplanar fluctuations, we preserve the lattice 
structure by coupling in a discrete fashion an infinite ar- 
ray of two-space and one-time QNLSM. We shall call 
this effective theory the quasi-two-dimensional QNLSM. 
The ratio c(3^p/^' is then determined by the parameters 
of the quasi-two-dimensional QNLSM, i.e., by the two- 
dimensional spin-wave velocity c and spin stiffness (or 
the gap A), as well as on J' /J. These are related to 
the microscopic parameters of the quantum Heisenberg 
AF. While the model with nearest-neighbor couplings 
and with physical spins (5 > 1/2) is known to be in the 
RC regime at low temperatures, the addition of terms in 
the Hamiltonian (|2.ip such as frustrating next-nearest- 
neighbor couplings, say, allows us to realize the QC or 
QD regimes. Wc thus consider the two-dimensional spin- 
wave velocity c and spin stiffness (or the gap A) as 
phenomenological parameters. 



The scaling functions we are looking for will depend 
on the ratios of the lengths cPaf, and, if we are inter- 
ested in an observable of the two-dimensional underlying 
system, of ^^^''{Paf)- However, at the transition tem- 
perature and in the limit of J'/J — > 0, only one of the 
two ratios is independent, and the scaling functions will 
indeed depend only on one dimensionless parameter 



B. Universal scaling functions 

The first universal scaling function, suggested in Ref.[^, 
involves the two-dimensional static and staggered spin 
susceptibility, 

^lin^^ J'xf^ iP^p) = F, (c/3^^/^(2) ^ (2.7a) 
where 

xi'^ (Paf) ■■= X^'' (kAF, ^ = 0; Paf) (2.7b) 
with kj^p ~ (tt, tt). In the RC regime, 

(cPaf/C^'HPaf)) = ^^i(O) ^ (2.8) 

MC simulations give C^^) ~ 1.3^ while the RPA approx- 
imation leads to C*-^-* = 2. Using the N = oo approxima- 
tion to the quasi-two-dimensional QNLSM, we will find 
(■(^^ = 1 in the RC regime. 

The two-dimensional spin-susceptibility is inaccessible 
experimentally at the three-dimensional ordering tem- 
perature. In order to mimic the two-dimensional AF 
wave vector in the three-dimensional system, we choose 
to work at the wave vector (tt, tt, 0). We claim the exis- 
tence of the universal scaling function 

^Ihn^^ J'xf^ {Paf) = ^2 {cPaf/S.') , (2.9a) 



where we have defined 

xi'' {Paf) ■■= X^'^ (tt, ^, 0, = 0; Paf) ■ 



(2.9b) 



Observe that any vector of the form (tt, tt, k^) with ^ 
TT would lead to the same conclusion but for different 
universal scaling functions that depend on k^ . In the N = 
oo approximation of the quasi-two-dimensional QNLSM, 
we find that (x) = 1/4 takes the same value in the RC, 
QC, and QD regimes. This suggests that this scaling 
relation is rather robust and well suited for numerical 
studies. 

Finally, we claim the existence of the universal scaling 
function 



(2.10a) 



■J'/ 



where S^'^^ is the instantaneous structure factor 

'•+°° d^im x(''(g,g.,^;/3) 



5(3)(q,^^,t = 0;/3) :=- 



(2.10b) 
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In the N = oo approximation of the quasi-two- 
dimensional QNLSM, we find that F^{x) ~ (a;/4) cothx. 

The reminder of the paper is devoted to proving the 
existence of these three universal scaling functions and to 
their computation in the N = oo approximation of the 
quasi-two-dimensional QNLSM. 



III. QUASI- TWO-DIMENSIONAL QNLSM 

The quantum nonlinear a model (QNLSM) was 
successfully used to study the low-energy and long- 
wavelength physics of the quantum one-dimensional 
Heisenberg AFji^ the quantum two-dimensional Heisen- 
berg model AFf^ii^ and the quantum quasi-two- 
dimensional Heisenberg AF.— In Haldane's mapping of 
the quantum Heisenberg AF to the QNLSMjH a crucial 
role is played by the correlation length. A large corre- 
lation length gives the possibility to separate slow from 
fast modes. Integration over fast modes can then be car- 
ried out perturbatively, retaining only slow modes. In 
the quasi-two-dimensional quantum Heisenberg AF, the 
length scale ^' defined in Eq. (|2.3ap can also be used as a 
characteristic length scale to separate the fast from the 
slow modes in Haldane's mapping. These fast modes can 
then be integrated out following the procedure used by 
Haldane,— and the partition function can be expressed 
in terms of a path integral over unit vectors. In the 
isotropic limit, the three-dimensional QNLSM, which is 
a pure field theory, is recovered. However, as noticed in 
Refs. 1113, S and O, the large values of allows us to 
take the continuum limit only within the planes. The 
resulting partition function has the form 



Z 



I I l[V[a^{x,r)]S{a'^{x,r)-l)l 



.-SM 



(3.1a) 

where (Tj{x,T) G M^, the integer = 3 for real spin 
systems, and r is the imaginary time. The action in 
Eq. (|3.1ap can be divided in two parts. 



S — Sq + Si, 



(3.1b) 



where 



J 



(3.1c) 

is the action of the two-dimensional QNLSM on a collec- 
tion of independent planes labeled by j and 



Si := 



2 J 



(3.1d) 



J 



describes the interplane couplingii^ The bare planar spin 
stiffness = JS^ and the bare planar spin-wave veloc- 
ity c = 2\/2SJ provided units such that a = 1 have been 



chosen. The 0{N) model is constructed allowing N to 
take any value larger than 2 in N. The local constraint 
cr'j{x,T) ~ 1 can be ensured by a Lagrange multiplier 
Xj{x,T). The N = oo approximation is then obtained 
after integrating out the original fields ctj and by ex- 
pressing the original partition function in the form 



(3.2a) 



where the effective action is now 



SA\]^^^^Y.jdrjd'^\ 



{x,t) 



J 



(3.2b) 



1 



Trln(-A,_,-aA, + A^-(a;,r)) 



2 ~' ^ —X,T Z ■ ■ 'J > 

We have introduced the bare coupling 

g cN/p^, 
the bare spatial anisotropy strength 

a:= J7J=(l/0^ 
(remember that a = 1), and 



(3.2c) 



(3.2d) 



(3.2e) 
(3.2f) 



The parameter N enters explicitly the action as a pref- 
actor only. In the limit of large- A^, any observable can 
be expanded in powers of 1/A with the leading order 
corresponding to the saddle-point approximation, [for a 
review, see Ref. [TH . The two-point spin correlation func- 
tion is 



Gf^,(a;,x') = 6''^Z-^ 
with 

G- ,(x,a;';A) {x,j 



I?[A]e-^^-«[^lG,,,(x,x';A), 



(3.3a) 



-A,., - «A, 



X6 



(3.3b) 



to leading order in the large-A^ expansion. We have in- 
troduced X = {x,ct). Because we approach the Neel or- 
dering temperature from above, we can assume isotropy 
in spin space of the spin correlation functions and drop 
the spin indices on both sides of Eq. p.3ap . The spin- 
isotropic susceptibility is then related to the correlation 
function in momentum spaced 

G{k,k,,uj) ^ / e-'(^-+'^-^^")G,-o(x,0) (3.4) 

■ J X 
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by 



'Nc 



G{k,k^,u). (3.5) 



In the saddle-point equation, which becomes exact in the 
hniit N oo, the ansatz = const is chosen in order 
to minimize the action S'gg[A]. This leads to the saddle- 
point equation 



1 



(3.6a) 



when 



.16 



gS5(x,x;AJ 



- y 



G(")(p,p„n7;AJ 



(3.6b) 



with 



p2 _|_ 2a{l — cosp^) + A^ 
(3.6c) 



G(°)(p,p„n7;A, 



The saddle-point equation p.6p is equivalent to the con- 
straint (cr|(x)) = 1. It can either be solved for any tem- 
perature that admits a nonvanishing and positive solu- 
tion A^ > 0, which is then related to the N = oo approx- 
imation of the temperature-dependent correlation length 



f = f (3) = J — 



(3.7) 



or with A^ = in order to determine the N = oo approx- 
imation of the critical temperature T^p = 1 //Sj^p ■ 

We will use the saddle-point equation (|3.6p to verify 
the existence of the three universal functions defined in 
Sec. Ill B] in the iV = oo approximation. 



NEEL TEMPERATURE IN THE TV = oo 
LIMIT 



IV. 



The saddle point (|3.6p is UV divergent. We follow 
Ref. H and use the relativistic Pauli-Villars regularization 
with cutoff A ^ c/3 for the propagator in Eq. p.6p . The 
form of the saddle point (|3.6p then depends on the ratio 
between the coupling g and its critical value in the two- 
dimensional QNLSM 



47r 

9c :=X. 



(4.1) 



When 1/g - 1/g^ > 0, the saddle point for the 

Neel temperature reduces to 



AF 



N 



(4.2a) 
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FIG. 1: (Color online) Plot of the auxiliary function Y{x) 
defined in Eq. (|4.4a|) as obtained by numerical integration. 



where we have introduced the renormalized spin stiffness 

Ps ■■= cN (---)■ (4.2b) 
\9 9c J 

If l/.g^ ~ 1/g > 0, the saddle-point ([SH) for the Neel 
temperature reads 



(4.3a) 



(4.3b) 



A /3ap = +2Y{cPaf/^') 
with the quantum disordered spin gap 



Attc I 

.9c 9 



in the two-dimensional QNLSM. In order to get Eq. (|4.2p 
and Eq. (j4.3p . we first performed the summation over 
Matsubara frequencies followed by the momentum inte- 
gration using the Pauli-Villars regularization scheme. In 
doing so we find the function 



Yix) 



dip 



In [2 sinh (a; sin ip^'\ 



(4.4a) 



which is monotonically increasing for .t > 0, has a first- 
order zero at .Tq ~ 0.93, and the asymptotes 



Y{x) ~ < 



Ina; 4- f^, a; < 1, 
yi{x — a;g), x — Xq, 
a; > 1, 



(4.4b) 



2x 7 



with the derivative yi = Y'{xq) ~ 1.22 and the constant 
7 ~ 0.53. We plot in Fig. [l]the function Y{x). As we 
shall show in the sequel, these three asymptotic regimes 
allow us to identify the three regimes (|2.4p . (|2.5p . and 
(|2.6p for which the Neel temperature is finite. 

The limit J' / J ^ always probes the RC regime (|2.4p 
provided the underlying two-dimensional subsystem is 
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FIG. 2: (Color online) Projection of the phase diagram of 
the quasi-two-dimensional QNLSM into the plane spanned by 
g/gc — 1 and J'/ J. The upper bold line is the phase boundary 
between the AF ordered and the paramagnetic phases corre- 
sponding to the fixed ratio R+ = n/4. The hatched regions 
denote the crossover between the QD and QC regimes, on the 
one hand, and between the QC and RC regimes, on the other 
hand. Each curve represents a possible trajectory with fixed 
R± used in order to construct the scaling functions in the 
limit J' / J 0; the temperature along the curves is always 
fixed to Taf- 



not quantum disordered, i.e., it satisfies the condition 
g < g^. The limit J' /J — > can only probe the QC 
regime (|2.5p or the QD regime (|2.6p after fine tuning. 
This is so because, for any value g/g^ ^ 1, the uncon- 
strained limit J' I J ^ brings the quasi- two-dimensional 
system either in the regime (|2.4p . that is dominated by 
two-dimensional classical fluctuations, when glg^<\ ox 
in the quasi-two-dimensional paramagnetic phase, that is 
dominated by two-dimensional quantum fluctuations at 
very low temperatures, when gj g^ > 1. When g/g^ < 1, 
the QC regime (j2.5p can only be probed in the J'/ J 
limit if the ratio of length scales 



R_ 



(4.5a) 



A. RC regime 

When g < gc and cf3^p ^ Eq. (|4.2[) simplifles to 



exp 



N 



(4.6) 



Observe that, for any flnite p^, the condition (|2.4p is al- 
ways met for J' small enough, as, in the limit J' / J 0, 
Eq. (|4.6p has the solution 



Paf{J'IJ) 



N 



2 I J' 



In 



cln(J/J') 



(4.7) 

with corrections of order jy jf) • The dependence of 
the Neel temperature on J'/ J shows an essential singu- 
larity at J'/J = 0; expansion (|4.7p is very poor. 



B. QC regime 



When g < g^ and c/3j^p/^' < Xq, Eq. (|4.2p becomes 



c(3 



AF 



cN 
47rp^e ■ 



(4.8) 



Equation (|4.8p is self-consistent with the assumption 
^PafI^' ^ ■^"O provided R_ ^ 1. The Neel tempera- 
ture can be expanded in a power series in 1/R_ when 
approaching the two-dimensional quantum critical point 
~ = from the two-dimensional quantum ordered side. 



as 



U g > g^ and cPj^p/S^' > Xq, Eq. (|4.3p can be rewritten 



cf3 



AF 



e 1 - l/i2y,R, 



(4.9) 



to be consistent with the assump- 
i,Q, we have to require R^ ^ 1. Again, 



In order for Eq. 
tion c/3j^p/C > Xq 
the Neel temperature can be expanded in a power series 
in 1/ when approaching the two-dimensional quantum 
critical point A = from the two-dimensional quantum 
disordered side. 



is held fixed. When g/g^ > 1, the QC regime (|2.5p or 
the QD regime (|2.6p can only be probed in the J'/J 
limit if the ratio of length scales 



c/A 



(4.5b) 



is held fixed. In the latter case, the magnitude of R_^_ 
distinguishes the fate of the limit J'/J 0, i.e., whether 
the onset of AF long-range order is characterized by the 
quantum fluctuations of regime (|2.5p or the quantum fluc- 
tuations of regime (|2.6[) . When R_^_ is below the thresh- 
old value of 7r/4, we shall see that three-dimensional AF 
long-range order is impossible down to and at zero tem- 
perature (see Fig. [2). 



QD regime 



When g > g^ and c/3^p 3> Eq. (|4.3[) simphfles to 



c/3 



AF 



27 



i?4 



Af ' 



(4.10) 



provided R^ > tt/4. When R^ < 7i'/4, the Neel tempera- 
ture is predicted to vanish in the N ~ oo approximation. 
The minimal interlayer coupling for which the three- 
dimensional Neel state exists can be obtained by solving 
Eq. ([iTTU)) at Pjip = oo. This gives 



J 



IT A 

17 



(4.11) 
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SCALING LAWS IN THE N ^ oo LIMIT 



FIG. 3: (Color online) Boundaries between the AF ordered 
and the paramagnetic phase as a function of the spatial 
anisotropy. The dependence on J' /J of the two boundaries 
(thick lines), one emanating from the origin {g,T) = (0,0) 
into the plane g — 0, the other emanating from the two- 
dimensional quantum critical point {g,T) — ((/c,0) into the 
plane T — 0, are given by Eqs. (|4.7|l and (|4.12|l . respectively. 
The curves on the critical surface (dashed lines) lead to the 
QC (lower line) and QD (upper line) behaviors. 



D. Critical coupling 



The critical coupling defined in Eq. (|4.2bp separates 
the Neel ordered phase from the paramagnetic phase of 
the two-dimensional system at zero temperature. We 
have seen that the quasi-two-dimensional ordered phase 
can exist in the limit J'/J — *■ even though g > g^, 
provided the two-dimensional quantum gap A is scaled 
accordingly. This implies that any finite coupling J'/J 
between planes renormalizes g^ to a larger value g^. The 
N ~ oo approximation for g^ can be obtained from 
Eq. (|4.10p after inserting in c/C — irA/A the renormal- 
ized gap A from Eq. (|4.3b[) . One finds the dependence 
on J'/J, 



1 

9c 



(4.12) 



We are now going to verify the claims (|2.7p . (|2.9p . and 
(|2.10p in the N ~ oo approximation. 



A. First scaling law 

One way to proceed in the case of Eq. (|2.7p is as fol- 
lows. We solve the saddle-point equation (|3.6p a first 
time with = but J'/J > to determine the Neel 
temperature, as was done in Sec. IIVI We then solve the 
saddle-point equation (j3.6p a second time but now with 
= J'/J = to extract the two-dimensional correla- 
tion length ^^^'(/3a_f) the Neel temperature, which in 

(2) 

turn determines Xs (Paf)- Alternatively, we can extract 
C(2)(/3^^) from 



of] I x, X 



, (5.1) 



where a matrix element of the TV = oo approximation to 
the propagator G'-^-' is explicitly given in Eq. (|3.6cp . This 
equation has the advantage of being cutoff independent 
in the Pauli-Villars regularization scheme. 

After summing over the frequencies and integrating 
over the two-dimensional momenta, Eq. (jS.ip reduces to 



= 2 



We insert this ratio into 



arcsinh ( -e^'^.'^^Ap/^') 
2 



(5.2) 



(5.3) 



where we made use of Eq. (j3.5p using the N ~ oo approx- 
imation of the Green function with a ~ 0. After trading 
J' for J and ^' defined in Eq. (|2.3ap . we find 



JXs [Paf) - [ ^, 



{cPafU'Y 



4 arcsinh ( e 



(5.4) 



(5.5) 



_____ 

Equation (|5.5p tells us that the function J'xs is a scaling 
fimction of the scaling variable cfiji^p/S^' . This is not F-^ 
from Eq. (|2.7p . since we are looking for a function of 
'^PafH^'^HPaf) ^^'^ its limiting value when J'/J —> 
has yet to be taken. 



Renormalized classical regime 



that is characterized by an essential singularity at J'/J = 
0. The boundaries (|4.7p and (|4.12p are depicted in Fig. [31 



In the RC regime (|2.4p . where c(3ji^p/SJ <C 1, and after 
expanding the right-hand side of Eq. (j5.2p in powers of 



8 



<^PafI^'^ we get 



and from Eq. (|5.5p 



24 V ^' 



(5.6) 



(5.7) 



According to Eq. ()5.6p . we can replace ^' by S.'^^HPaf) 
Eq. (|5.7|) to the first nontrivial order, 



AF 



12 K^^KPaf) 



(5.8) 



(2) 

We conclude that, in the N = oo limit, J'xs {(^af) 
equals the scaling function 



with the scaling variable 



c(3 



AF 



e'HPAF) 



(5.9a) 



(5.9b) 



in the RC regime p.4p . The universal scaling function Fi 
is then obtained from JF]^^(x) taking the limit J' /J 0. 
Observe that we do not expect the function J^^'~^{x) to 
be universal in general. For example, this function is 
modified by adding longer range interlayer couplings in 
Eq. Odl) ^ 

For comparison with numerical simulations, it is in- 
structive to compute this relation as a function of a = 
J'/ J. Inserting the Neel temperature from Eq. (|4.7|1 in 
Eq. (jO)) leads to 



J'xi^)(/3^^)^1- 



1 C^N^ J' , 2 f J 



■In" 



12(47rpj'J V^' 



(5.10) 



The universal value in the limit J'/,/ ^ is thus estab- 
lished, since 



F,{Q) = 1 



(5.11a) 



is independent of any microscopic details, i.e., the spin 
stiffness or the spin-wave velocity c. For comparison, 
the RPA predicts 



F^^^^\0) = 1/2, 



(5.11b) 



while the MC calculations from Ref. gives 

Fi^^^^(O) -0.77. (5.11c) 

We show in Appendix |A] that the form of the nonuni- 
versal corrections in Eq. (j5.10|) are similar to the correc- 
tions obtained in the quasi-one-dimensional Ising model. 
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FIG. 4: (Color online) The scaling function Fi is plotted vs 
J' /J for spin 1/2. The curve represents the N — oo ap- 
proximation (|5.8|l . with the Neel temperature obtained from 
numerical solution of Eq. ([4.6|l . The points are the MC data 
from Yasuda et al. in Ref. In both cases, the normalization 
is such that the curves go through 1 for J'/ J = 10~^. 



The factor cN/{Attp^) scales as \/S. As expected, the 
nonuniversal corrections get smaller for bigger spins. We 
show in Sec. IVII that the nonuniversal corrections vanish 
in the N = oo approximation to the classical Heisenberg 
AF. 

The result (|5.10p is only reliable for very small J'/ J, 
due to the limited quality of the approximation (j4.7p . 
To compare the N = oo approximation with MC simula- 
tions, it is best to rely on a numerical solution of Eq. (14. 6p 
for the Neel temperature. As can be seen from Fig.|4l for 
J' /J < 0.1, the N = oo approximation (|5.8p of the first 
scaling law obeys a limiting behavior as J' / J sim- 
ilar to the MC results from Ref. 0. For J'/ J > 0.1, the 
magnitude of the deviation of the effective coordination 
number from its universal asymptotic value as J'/ J ^ 
can be explained using the empirical formula for the Neel 
temperature proposed in Ref. |j. However, the sign of the 
nonuniversal correction in Eq. (|5.8p is wrong. We do not 
expect the quasi-two-dimensional QNLSM used here to 
be valid for the description of the system at J'/ J > 0.1 as 
the effective interlayer separation ^' is only a few lattice 
spacing large. Instead, it is likely that it is necessary to 
also consider the lattice structure of the system in each 
plane to obtain an appropriate description. 



2. Quantum critical regime 

The condition A = fine tunes the underlying two- 
dimensional system to be at a quantum critical point. 
From Eqs. ((13)) and (j4.4bp . we then have cP^p/^' = x^. 
Taking advantage of this relation in Eq. (|5.2p gives 



e'Kf3AF) 



e_ 



(5.12a) 
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where 



4- Summary 



l + \/5 



e := 21n 



The right-hand side of Eq. (|5.4p thus becomes 



02- 



Chubukov et al. in Ref. |8| have shown that 

cl3 



= e 



(5.12b) 



(5.13) 



(5.14) 



in the N = oo approximation when A = 0. The universal 
function from Eq. (|2.7p thus takes the value 



(5.15) 



exactly at the QC point A = 0. 

In the scaling limit J'/ J ^ with the ratios from 
Eqs. (|i3a|) and ((45bl) held fixed, Eqs. (jSl^ and (jSl^ 
become 



and 



respectively, where 



(5.16) 



(5.17) 



(5.18) 



3. Quantum disordered regime 



We close with the QD regime (|2.6p for which 
cf3j^p/^' 3> 1. The underlying two-dimensional system 
is in its quantum disordered phase with 



C^^^T = 0) = c/A. 



(5.19) 



Taking the scaling limit J' /J with the ratio from 
Eq. dfSbJ, 

= ^^'HPAF)/e + 0[e-^^'--/^), (5.20) 
held fixed turns Eq. (|5.4p into 

(5 21) 

up to exponentially small corrections in (R, — n/A)^'^/^. 



To sum up, we have the N = oo approximation 
' 1, RC regime. 



^nm^J'x(''(/3^j.) 



QC regime. 



QD regime, 

(5.22) 

to the first universal function (|2.7p . where R-^ ^ 1 in 
the QC regime while R^ ^ 7r/4 in the QD regime. The 
result F^(0) = 1 is exact in the RG regime !^7r\ . The 

result {cP af I ^^"^^ af)) = q2 — is accurate up to 
corrections of order in the QC regime (|2.4p . The 
result F-^{cl3j^p/x^'^\l3j^p)) ~ R\ is accurate up to ex- 
ponentially small corrections in (i?_|_ ^ 7r/4)'^/^ in the 
QD regime ()2.5p . The RC regime ()2.6p is the only one 
for which we were able to compute the nonuniversal cor- 
rections to the limit J'/ J 0. They are of second order 
in the variable cPj^p/^^'^\(3j^p) but nonanalytic in the 
variable J'/ J, namely of order {J' / J)\r?{J / J'). 



B. Second scaling law 

The second universal scaling law (|2.9p . which involves 
the three-dimensional spin susceptibility, is established 
as follows. In the N = oo approximation. 



S'g 1 



(5.23) 



cN 4a' 

which immediately leads to the second universal relation 

J'xf\f3AF)^\, (5.24) 
in the RC, QC, and QD regimes. 



C. Third scaling law 

We close by establishing the third universal scaling 
law (|2.10p . which involves the instantaneous structure 
factor. We insert the definition (|3.5p in Eq. (j2.10bp . and 
use the N = oo approximation to the correlation func- 
tion. After performing the frequency integration, we find 

Sik^p, k, = 0, i = 0; 13) = l^e' coth {cf3/e) ■ (5.25) 
The last universal relation is then established, 

^ApJ'SikAF. k^=0,t = 0; ^^p) = ^3 (c/3^^/e') , 

(5.26a) 

where in the N = oo approximation. 



J^^{x) := -cothx 



(5.26b) 
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and 



0(3 



AF 



(5.26c) 



The universal scaling function F3 is then obtained from 
!F^{x) taking the limit J'/ J ^ 0. As for the first scaling 
law, we do not expect the function J-^ {x) to be universal 
in general. The behavior of J-^{x) when x <C 1, x ^ 1, 
and a; 1 distinguishes the RC, QC, and QD regimes. 
Observe in this context that, as T^{x) — > 1/4 + 0{x^) 
in the limit of small x, the function i^g is the same as 
F2 in the RC regime, as expected from the quasielastic 
approxiniatio] 



.17.18 



x(g,c^ = 0;/3)~-/35(9,t = 0;/3). 



(5.27) 



In the QD regime, the scaling variable x can become ar- 
bitrary large as i?+ 7r/4. [In particular, this is the case 
when J' J^, see Eq. (|4.1ip .] Consequently, the func- 
tion .^3(0;) ~ x/4 on the right hand side of Eq. (|5.26ap 
diverges when the ratio that defines how the limit 
J' ^ is taken approaches 7r/4. 



VI. SCALING LAWS IN THE N = S = 00 LIMIT 

The three components of the spin operators with spin 
quantum number S commute in the limit S* — > cxd. The 
results of Yasuda et al. suggest that the effective coor- 
dination number (|1.2p is independent of S in the limit 
J' / J 0. This limiting value of the effective coordi- 
nation number is thus the same irrespective of S taking 
values deep in the quantum regime, S = 1/2, or clas- 
sical regime, S = 00. If so, it is instructive to compare 
our analysis of the quasi-two-dimensional QNLSM (j3.1bp 
with that of the quasi-two-dimensional NLSM. The ac- 
tion of the quasi-two-dimensional NLSM is obtained from 
Eq. (j3.1b[) by dropping all references to the imaginary 
time, setting S — 1 and c = 1, and replacing the integral 
over T by 1/T, 



2T ^ 

j 



J' 

T 



(6.1) 



with = J. Again, the constraint cr']{x) ~ 1 can be 
removed in favor of the Lagrange multiplier A. If the 
original fields cr are integrated out, one obtains the effec- 
tive action 



x\j + ]^ Tr In ( - - a + A^- ) 



where 



(6.2a) 



g = NIp^ (6.2b) 

that should be compared with Eqs. (|3.2b|) and (|3.2cp . 
The spin susceptibility is now related to the two-point 
correlation functions by 

x('Hfc + fcAF,7r + fe;/3) = j^G{k,k,;P) (6.3) 



that should be compared with Eq. (|3.5p . 

To compute the first scaling law in the N = 00 ap- 
proximation, we proceed as in the quantum case. We 
thus solve directly the classical counterpart to Eq. (|5.1|) 
for the ratio C'/^'^' '^ith ^^^^ evaluated at the Neel tem- 
perature. We find 



e^KPAF) 



1. 



(6.4) 



Hence the first scaling law in the N ^ 00 approximation 
simply reads 



with 



and 



J'xf\pAF)^^li^) 



^i(x) = l 



JPafO- 



(6.5a) 



(6.5b) 



(6.5c) 
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The counterpart to the universal function (|2.7p obtained 
in the limit J'/J ^ is then 



F,ix)^l. 



(6.6) 



It agrees with the RC limit (|5.11ap . This is consistent 
with the MC study from Yasuda et al. in Ref. [3. Con- 
trary to the quantum case, nonuniversal corrections are 
absent in the saddle-point approximation. This result 
is also consistent with the observation from Ref. that 
the finite J'/J nonuniversal corrections decrease with in- 
creasing 5*. 

We turn to the second scaling law. Inserting the = 
00 approximation of the two-point correlation function 
in Eq. (|6.3p at the wave vector (tt, tt, 0) leads to 



AF) 



(6.7) 



as in the quantum case. 

There is no independent third scaling law in classical 
thermodynamic equilibrium since all obscrvablcs are time 
independent by assumption. 



VII. COMPARISON OF THE SCALING LAWS 
WITH EXPERIMENTS 

Many quasi-two-dimensional quantum Hcisenberg AF 
are now available. Neutron studies have been performed 
on compounds with spin S = 1/2 [La2Cu04 in Refs. [19I 
and [10, Sr2Cu02Cl2 in Refs. [21] and [2^. and copper 
formate tetradeuterate (CFTD) in Ref. |23| spin 5* = 1 
(La2Ni04 in Ref. HJ and K2NiF4 in RefHH), as well as 
spin S = 5/2 (KFeF4 in Ref.lHand Rb2MnF4 in Refs.lH 
and [13). All these materials are good realizations of the 
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two-dimensional QNLSM above their ordering tempera- 
ture Tj^p. In particular, the temperature dependence of 
the correlation length in these systems is well explained 
(without free parameter) by the correlation length of the 
two-dimensional QNLSM in the RC regime and at the 
three-loop order fit^.!^^ 



c/a 
211 ps 



,27rp,/T 



1 - 



1 



T 



(7.1) 



All these systems have an interplane exchange coupling 
J'/ J which is much smaller than the spin anisotropy. 
The latter is either of XY type (for the 5 = 1/2 systems 
listed above) or Ising type (for the remaining examples 
chosen) . Hence the onset of AF magnetic long-range or- 
der at Tj^p is a classical critical point that belongs to 
the XY or Ising universality class. Moreover, T^^, is 
mainly determined by the (lower) critical temperature 
of the two-dimensional system, which is finite in opposi- 
tion to that for the spin-isotropic two-dimensional quan- 
tum Heisenberg AF. For this reason, we do not believe 
that the scaling laws derived above for a spin-isotropic 
quantum Heisenberg AF can be observed in the materi- 
als mentioned above. 

A family of organic quasi-two-dimensional quantum 
AF has been synthesized with the general chemical for- 
mula A2CUX4, where A=5CAP (5CAP stands for 2- 
amino-5-chloropyridinium) or A=5MAP (5MAP stands 
for 2-amino-5-methylpyridinium), and X=Br or Cl^ 
The AF Heisenberg exchange coupling J in these ma- 
terials is between 6 and lOK, and they have a Neel tem- 
perature of about Taf — 0.6J. According to Ref. H^, 
the dominant sublcading term in the Hamiltonian de- 
scribing their magnetic response is the interlayer Heisen- 
berg exchange coupling J'. The rationale for J' being 
the leading subdominant term to the planar J is the fol- 
lowing mean- field argument 3^ Consider the compound 
La2Cu04 as an example of a layered structure in which 
adjacent layers are staggered, i.e., there exists a relative 
in-plane displacement (1/2, 1/2) between any two neigh- 
boring layers. Each in-plane ion has thus four equidistant 
neighbors in the layer directly above it. Assume that, in 
any given layer, the spin degrees of freedom occupy the 
sites of a square lattice and that they are frozen in a clas- 
sical Neel configuration. Any in-plane spin has then four 
nearest-neighbors in the layer above it whose net mean 
field vanishes. In contrast, the stacking of planes is not 
staggered in the organic quasi-two-dimensional AF de- 
scribed in Ref. [2Sl . so that this cancellation mechanism 
does not work. Consequently, the ratio J'/ J is expected 
to be much larger for A2CUX4 than for La2Cu04. 

The interlayer exchange parameter for A2CUX4 is esti- 
mated in Ref. [H to be J' /J ~ 0.08 with the help of the 
formula 



(7.2) 



Equation (|7.2p can be found in Ref. while its quasi-one- 
dimcnsional version was obtained by Villain and Loveluck 



in Ref. 'so'. The Neel temperature of the spin-1/2 quasi- 
two-dimensional Heisenberg model is hereby estimated 
by balancing the thermal energy T^p against the gain 
in Zeeman energy obtained by aligning a planar spin 
S{S + 1) along a mean-field magnetic field of magnitude 
zJ'[£^^^'>{T^p)/a]'^ that points parallel to the stacking di- 
rection. Here ^'•^•'(Taf) is the two-dimensional correla- 
tion length estimated at the Neel temperature from for- 
mula (|7.1|) and z is the coordination number of layers. 
Using their empirical formula (which is consistent with 
the first scaling law in the limit J' /J — > 0), Yasuda et al. 
in Ref. claim the three times larger value J'/ J ~ 0.24. 
This discrepancy is explained in Appendix [Bl To our 
knowledge, there is no independent estimate for the in- 
terlayer coupling J'/ J deduced from measuring the spin- 
wave velocity for A2CUX4. We are also not aware of 
measurements of the correlation length above the critical 
temperature for A2CUX4. 



A. Experimental discussion of the first scaling law 

The two-dimensional spin susceptibility is inaccessible 
to a direct measurement at the Neel temperature, the 
temperature at which our scaling laws hold. This is not 
to say that the two-dimensional spin susceptibility is in- 
accessible at all temperatures. In fact, one would expect 
a window of temperature above the Neel temperature for 
which the three-dimensional spin susceptibility is well ap- 
proximated by the two-dimensional one in a good quasi- 
two-dimensional AF magnet. If so one could try to test 
the first scaling law by extrapolation from high temper- 
atures. 

To this end, we could first attempt to use the measured 
values of the Neel temperature and of the anisotropic 
spin-wave dispersion together with computations of the 
two-dimensional spin-wave velocity c and spin stiffness 

from first principle, to test the accuracy of the predic- 
tion (15.611 for the ratio 



(7.3) 



between the two-dimensional correlation length 
^'■^'(r^^) and the effective interlayer spacing (|2.3ap . 
Indeed, it is reassuring to know that there exists a good 
agreement between formula (|7.ip and measurements of 
the correlation length above Tjnp in the RC regime. 

Needed are the numerator and denominator in the RC 
regime of the right-hand side of Eq. (|7.3[) expressed as a 
function of c, p^, and T^^, to first order in the 1/A^ 
expansion. The two-dimensional correlation length is^ 



^^^>{T) _ c/a ((N-2)T 



l/{N-2) 



exp 



27rp, 



(iv- 2)r' 

(7.4a) 



where the constant 

eo = (e/8)'/(^'''r(l + l/(A^-2)) (7.4b) 
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depends on e = 2.718... and the Gamma function. The 
Neel temperature is the imphcit solution to^ 

.ln4fe + ^ln. '^'^ 



{N-2)T^p cya^ N-2 {N~2)T^p 

(7.5a) 



1.012 1 J 
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N-2 Z' ,/' 



with 



37r2iV IQpJ 



X 1 + 



1.069\ ({N -2)T 



N 



l/(N-2) 



.(7.5b) 



The dependence on the momentum cutoff A of the two- 
dimensional correlation length occurs through ^' only, 

One thus gets the universal number 

,1.506 \ 1/(^-2) 



s = r(i + i/(iv-2)) 



(7.7) 



as it is independent of the spin-wave velocity and the 
spin stiffness. In particular, S « 1.127 for N = 3. (In 
the = oo, we found S = 1.) 

We have computed the ratio S for some spin S = 
1/2 quasi-two-dimensional AF using the values of T^^p, 
(J'/J)s^, 27rpg, and c hsted in Table HI Here, {J'/J)sw 
is the ratio of the interplane to the intraplanar exchange 
couplings deduced from the spin-wave spectra measured 
using inelastic neutrons scattering at temperatures well 
below the measured T^p. The measured value {J' / J)sw 
is interpreted as the multiplicative renormalization 



8 cAAf\ 
3^ 16^/ 



(7.8) 



that arises solely from quantum fluctuations,^" In turn, 
this allows us to express ^' in Eq. (|2.3ap in terms of the 
microscopic parameters a, A, J, and J' on the one hand, 
and the macroscopic parameters and Tj^p on the other 
hand, 



a 



'^'^ ^ 1.069^ ^''^ 



{N - 2)T 
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A^ 

-l/[2(JV-2)] 



(7.9) 



The crossover length scale ^' in Table U then follows 
from inserting A^ = 3 and the values for T^p and Pg 
from Table HI The same is done with Eq. (|7.4ap to 
obtain ^^^H^af) in Table [H One finds the measured 



values ^„ 



0.49 for La2Cu04, S„ 



0.02 for 



SraCuOaCla, and S„^^, w 0.34 for CFTD. The small- 
ness of the measured Sjjjcas compared to S « 1.127 in- 
dicates that the actual transition to the ordered phase 
takes place when the two-dimensional correlation length 
is smaller than the effective interlayer separation This 
discrepancy with the large-A'^ expansion can be under- 
stood as follows. The compounds of Table U all have 
a smaller anisotropy {J' /J)sw compared to the XY 
anisotropy or to the Dzyaloshinsky-Moria term. Further- 
more, the Sr2Cu02Cl2System has about the same critical 
temperature as the other two compounds, although its 
anisotropy is three orders of magnitude smaller. This in- 
dicate a^^i20,22,23 ^j^g^^ ^j^g pi^^gg transition is not triggered 

by a pure dimensional crossover but by a combination of a 
symmetry and dimensional crossover that effectively en- 
hances the true two-dimensional correlation length over 
that for the pure two-dimensional 0(3) QNLSM. The 
same conclusions may be drawn for systems with higher 
spins. 

We close this section by evaluating the first scaling law 
in the RC regime to the order 1/A''. The two-dimensional 
spin susceptibility at the wave vector (tt, tt) and at van- 
ishing frequency = is 



(7.10a) 



where the two-dimensional multiplicative renormaliza- 
tion Z is§^ 



(7.10b) 



X 1 + 



0.188\ f{N~2)T^p 



l/{N-2) 

N J \ 2np, ) 
Remarkably, the momentum cutoff A drops from the ratio 

2 I I 0.188 



N 



Z' 1 



1.069 

N 



X 2-1/(^-2). 



(7.11) 



For A^ = 3, Z/Z' = 0.391. (In the A^ = oo limit, we find 
Z / Z' = 1.) The 1/A'' correction to the first scaling law 
thus reads 



J'x^Ht^f)-^;^'- 



(7.12) 



For A^ 3, J'xs^\Taf) = 0.497. Observe that this is 
very close to the RPA prediction of 1/2. 



B. Experimental discussion of the third scaling law 

What has been extensively studied experimentally is 
an inelastic neutron-scattering measurement that is be- 
lieved to yield an approximation to the instantaneous 
structure factor 



-t-oo 



duo' 
"2^ 



S'o(7r, TT, 0,a;';/3). 



(7.13) 
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La2Cu04 


Sr2Cu02Cl2 


CFTD 




Refs. 19 and 20 


Ref. 22 


Ref. 23 


J 


15667^ 


1450-?^ 


73.27^ 




2074 7 


1769 7 


0.225J 


'J XY 1 'J 


5.7 X 10~* 


5.3 X lO"** 




.lr~iAT I.I 
'J U 1\1 1 " 


1.5 X 10~^ 




7 X 10~^ 


(J' /J) 


5 X 10~* 


10"* 


5 X 10~^ 


2tvPs 


1.131J 


1.15J 


1.31J 


c 


1.669Ja 


1. 669 Ja 


1.669Ja 




105.6a 


303a 


69.77a 




215.81a 


1.548 X lO^a 


207.197a 



TABLE I: Relevant parameters for some 5 = 1/2 quasi- 
two-dimensional AF. The intraplanar nearest-neighbor AF 
Heisenberg exchange coupling is J, while J' is the AF Heisen- 
berg coupling between nearest-neighbor layers. The XY ex- 
change coupling is Jxy, while Jum is the coupling of a 
Dzyaloshinsky-Moriya interaction. The Neel temperature is 
Taf- The two-dimensional spin stiffness 27rpa and spin- wave 
velocity c are obtained either from large S expansions or 
from MC simulations. The two-dimensional correlation length 
C^'^H^af) is obtained from formula (|7.1|l . Lengths are mea- 
sured in units of the lattice spacing a. 



Two ranges of temperatures have been studied in the 
literature. In the 1970s the chosen temperature range 
was a very narrow one about the AF transition tempera- 
ture. The rationale for this choice was to study the criti- 
cal regime surrounding the ordering temperature. In the 
1990s the temperature range was broader and above the 
onset of three-dimensional critical fluctuations fi^~— The 
rationale was primarily to study two-dimensional fluctu- 
ations associated with the classical renormalized regime 
of a two-dimensional AF. In either cases, the experimen- 
tal measurement of Sq{P) is performed in arbitrary units, 
i.e., the overall scale of 5q(/?) is unknown. Since we are 
predicting a universal number when measured in some 
given units, we need to convert any measured number 
in arbitrary units SQ^'"\f3) into a number in some cho- 
sen units. We do this by multiplying Sq^'^\/3) with a 
conversion factor, 



-,(a.u) 



(/3)x 



S, 



(MC) 



(7.14) 



obtained by taking the ratio of the number S^'^\f3*) 
at inverse temperature /3* computed within some scheme 
with the measured number in the same arbitrary units 
Sq'^'^\P*) at the inverse temperature (3*. 

We arc only aware of thepublished MC computation 
by Kim and Troyer in Ref. [32 of the static structure fac- 
tor at (tt, tt) of a spin-1/2 AF on a square lattice. We 
use this computation at temperatures such that the cor- 
responding correlation length is 3 in units of the lat- 
tice spacing. These are high temperatures for which 
the two-dimensional static structure factor should be a 



good approximation to the three-dimensional one. We 
are not aware of a published calculation of the static 
structure factor using either MC or high-temperature 
series expansion for spin-1 or spin-5/2 AF on a square 
lattice that are needed to reinterpret from our point of 
view the experiments from Refs. [13 or Refs. [1^, [26l . and 
[27l . respectively. If we restrict ourselves to the experi- 
ments on quasi-two-dimensional spin-1/2 from Refs. [lol 
[22[ . and[23l. we fail again to observe a signature of univer- 
sality in J'(3Sq{(3) when measured at the temperatures 
T ='337/C, 278/^, ISK that are 12K, 21.5K, and 1.5K 
above the corresponding ordering temperatures, respec- 
tively, as we find that J'(3Sq{(3) takes the values 1.17, 
0.001, and 0.077, respectively. 



VIII. CONCLUSIONS 

As we have seen, interesting scaling relations can be es- 
tablished between the interplane coupling J' and observ- 
ables of the underlying two-dimensional system or the 
quasi-two-dimensional system. The large- A'^ approxima- 
tion is particularly well suited to compute the universal 
functions obtained in the limit J' ^ at the Neel tem- 
perature and derive the leading nonuniversal corrections. 
The saddle-point approximation already leads to predic- 
tion in qualitative agreement with Monte Carlo (MC) 
simulations. To compare with existing experiments, we 
have included corrections of order 1/A^. These correc- 
tions allow us to take into account the renormalization 
of the wave function as well as the renormalization of 
the interlayer coupling, but will probably not affect the 
qualitative picture obtained in the saddle-point approxi- 
mation. 

The analysis of the quasi-two-dimensional model has 
revealed the existence of different regimes at low tem- 
peratures, corresponding to the renormalized classical 
(RC), quantum critical (QC), and quantum disordered 
(QD) regimes of the two-dimensional underlying system. 
The universal constant obtained from the first and third 
scaling laws in the limit J'/J are strongly modi- 
fied depending on the regime considered. For the second 
scaling law, we did not find any distinction between the 
three regimes in the saddle-point approximation. How- 
ever, we expect different corrections for the RC, QC, and 
QD regimes to the next order in the 1/iV expansion. 
While the RC result will probably not be strongly af- 
fected by renormalizations, the constant obtained in the 
quantum regimes might be modified more consequently 
by higher orders of the 1/N expansion. 

The first scaling law in the RC qualitatively agrees 
with the MC simulations of Yasuda et al. in Ref. [j for 
small J'/J < 0.1. Small variations of the Neel temper- 
ature are exponentiated in the correlation length. As 
a consequence, the first scaling law is very sensitive to 
small errors in the expression for the Neel temperature, 
so that a quantitative agreement between numerical sim- 
ulations and our results seems difficult to achieve, even 
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going beyond the saddle-point approximation. 

Recent MC simulations of quasi-two-dimensional sys- 
tems in the RC, QC, and QD regimes have been per- 
formed by Yao and Sandvik in order to test the scaling 
relations described in the present paper Their numer- 
ical results share many qualitative properties with our 
predictions. 

It is instructive to view our first scaling law against 
the mean-field estimate for the Neel temperature^i^S ob- 
tained by balancing the gain in energy derived from 
aligning planar spin S with a mean-field magnetic field 
S'J'[^(^)(T)]^ parallel to the stacking direction with the 
thermal energy T. This would give the naive estimate 



J 



rpMF 

J ■ 



(8.1) 



Instead, we can improve this estimate by including two 
types of logarithmic corrections. We can replace the bare 
ratio J'/J by (J'/J)g„ and we can multiply the right- 
hand side of Eq. (|8.ip by logarithmic corrections in the 
two-dimensional correlation length C}-"^^ (T) measured in 
units of the lattice spacing a. The Neel temperature is 
then the solution to 



T 



AF 



J 



ln« 



^^""{T^f) 



with the exponent q measuring the strength of the log- 
arithmic corrections. Note that while Eq. (|8.2I gives the 
correct scaling, it may have multiplicative nonuniversal 
corrections as exemplified by the presence of the lat- 
tice constant a on the left-hand side. For the case of 
the 0{N) symmetry classes with N > 2, wc have in 
fact shown the existence of logarithmic corrections with 
a quantum origin that are induced by the substitution 
J'/J {J' I J)svr of logarithmic corrections with a 
classical origin through the validity of Eq. (|8.2p with the 
extrapolation 



N - 1 



9 = 



N 



(8.3) 



of our leading (q = 1) and first subleading calculation 
((7 = 1 + 2/iV) in the 1/A^ expansion. It would be inter- 
esting to investigate the symmetry classes N = \ (Ising) 
and N = 2 (XY) from this point of view. 

The most important limitation to experimental ob- 
servations of the scaling laws proposed here are sym- 
metry crossovers. Indeed, most of the candidates 
for quasi-two-dimensional antiferromagnets have a spin 
anisotropy which is larger than the anisotropy in space. 
The anisotropy in the spin space leads to a crossover 
from an 0(3) Heisenberg model to an Ising or XY 
Heisenberg model which dominates over the dimensional 
crossover. The organic compounds (5MAP)2CuBr4 and 
(5CAP)2CuBr4 seem to be promising candidates for a 
magnet in which the dimensional crossover dominates 
over the symmetry crossover It would therefore be in- 
teresting to have independent estimates of the interlaycr 



coupling constant J'/J from the spin- wave dispersion on 
the one hand as well as from a measurement of the two- 
dimensional correlation length on the other hand for this 
class of compounds. 



Acknowledgments 

We would like to thank Daoxin Yao and Anders Sand- 
vik for showing us their numerical simulations before 
publication and Hcnrik R0nnow for useful discussions. 
M.B.H. was supported by U. S. DOE contract W-7405- 
ENG-36. 



APPENDIX A: QUASI-ONE-DIMENSIONAL 
ISING MODEL 

The first universal relation proposed in Refs. and 
can be established exactly for the strongly anisotropic 
two-dimensional ferromagnetic Ising model using the On- 
sagcr solution (sec, for example, Rcf . [ssh . 

Similarly to the two-dimensional Heisenberg AF, the 
one-dimensional Ising chain does not order at any finite 
temperature. Using the matrix formalism, the spin sus- 
ceptibility of the Ising chain can be computed exactly, 



+ 1 



tanh( J/?) ei/C - 1 ' 
where the correlation length is 

1 



lntanh(J/3) ' 



(Ala) 



(Alb) 



The two-dimensional Ising model has a nonvanishing 
transition temperature, which can be obtained as the so- 
lution to 



sinh (2/3c J) sinh (2/3^7') = 1, 



(A2a) 



where J and J' are the nearest-neighbor exchange cou- 
plings along and between the Ising chains, respectively. 
The solution of this equation is approximately 



— In In In — . 

2J J' 2J J' 



/3c w In ^ - ^ In In . (A2b) 
Using these results, one can compute the first scaling law, 

J'xi'HPc)=^iiJPc), (A3) 

where the scaling function J-i can be evaluated exactly. 
In the limit of small J', we obtain 

J^i (x) « 1 + 2e-2^, x = Jl3^. (A4) 

As a function of J' / J , the scaling relation reads 



1 J\ 2 J 

1 H In — . 

J J' 



(A5) 



This gives a value Q ~ \, which is half the coordination 
number expected in RPA. 
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APPENDIX B: 



TWO DIFFERENT ESTIMATES 

FOR Taf 



The interlayer coupling J'/ J can be estimated from 
the Necl temperature, provided the dependence of T^^p 
on J'/ J is known. We will discuss here different approx- 
imations giving this dependence. 

In the RPA from Rcfs.[3 and [3, the effect of the inter- 
layer coupling is encoded by an effective magnetic field 
proportional to the staggered magnetization in the ad- 
jacent layers that is multiplied by J' . The value of the 
staggered magnetization is determined self-consistently. 
The RPA staggered spin susceptibility is found to be 



^(RPA)(-y^ _ 



xf\T) 



1-zJ'xf (T) 



(Bl) 



where z is the coordination number of a single layer. At 
the Neel temperature, the spin susceptibility diverges. 
This condition leads to the following equation for the 
Neel temperature in the RPA 



AF 



1. 



(B2) 



This equation was the motivation for the first scaling 
law. It turns out that the quality of this approximation 
is ensured by the first scaling law. 

Another estimate for the Neel temperature in the 
quasi-two-dimensional Hcisenberg AF can be obtained by 
comparing the thermal energy Tj^p with the interaction 
energy between ordered spins in adjacent layers This 
argument leads to Eq. (|7.2p . 



zJ' 



M 



T 



AF 



(B3) 



where M/Mq ~ 0.36 is the reduction in the staggered 
magnetization at T = due to the two-dimensional spin 
fluctuations at length scales shorter than ^^^^ (T) . A sim- 
ilar equation has been derived for quasi-one-dimensional 
systems by Villain and Loveluck in Ref. [13 based on a 
real-space decimation argument. Equation (jBSp is in 
contradiction with the first scaling law that leads to 
Eqs. (ISH) and (HH). 

An estimate of the interlayer coupling J' / J from the 
Neel temperature and from the two-dimensional correla- 
tion length extrapolated down to Tj^p can be obtained 
using Eq. ((7?^ . Inserting ^' from Eq. ([7^ in Eq. (fTB]) 



leads to 



K 



^^'Ht^f) 



l/(W-2) 



XN-2)T^p^ 

(B4) 

where K = ^^/(l + ^^^)- For instance, using 2ttp^ = 
1.131J, C'-^HT^f) ^ 3.3a, and K = 0.937 for 3 leads 
to 

(J7J),„~0.32 ^ (B5) 
for the organic compounds A2CUX4 with T^p = 0.6J. 
This result is bigger than the estimate J'/ J ~ 0.24 pro- 
posed by Yasuda et al. in Ref. using their phenomeno- 
logical formula based on MC simulations, which itself is 
larger than the estimate J'/ J ~ 0.08 from Woodward et 
al. in Ref. [l^ based on Eq. (|7.2p . However, in view of 
the relative large J' / we expect large corrections to our 
scaling laws. 

As compared to the mean- field result (|B3[) . Eq. (|B4p 
contains logarithmic corrections. Indeed, using 



In 



(A^-2)T^^ 
we can rewrite Eq. (|B4p as 



O 



In In 



(2) 



T 



AF 



In"? 



(B6) 



where 



q 



N -I 
N -2' 



(B7) 



(B8) 



In particular, g = 1 in the N = 00 approximation, while 
for = 3 we find q = 2. The same value q = 2 can be 
obtained in the classical case, using the results of Brezin 
and Zinn-Justiui^ Indeed, they found in a 2+e expansion 
up to two loops, 



xl2)(T)~ r3gxp(47r J/T), 
for the staggered spin susceptibility and 

^(2)(T)^rcxp(2^j/r) 



(B9a) 



(B9b) 



for the correlation length. Equation (jB7p with g = 2 then 
follows from the first scaling law (|2.7ap . 
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